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Abstract. We obtain the viscosity solution of G-heat equation with the initial 
condition 4>(x) = x n for each integer n > 1 using the method of G-Brownian 
motion. 
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1 Introduction 

The notions of G-normal distribution, G-expectation and G-Brownian motion 
were firstly introduced by Peng (see [10] and [IT]) via the following G-heat 
equation: 

^~ G fe) = ' «lt=o = 6 (f,aO€[0,oo)xR, (1) 

where G(a) = \{a + — a 2 a~). Here a 6 [0,1] is a fixed constant and a + = 
max{0, a}, a~ = (— a) + . This equation is also the Barenblatt equation except 
the case a = (see Q], [5] and [5]). 

Under the sublinear framework, Peng gave many important notions corre- 
sponding to linear case and obtained many important properties of G-Brownian 
motion (see [13]). Recently, Peng developed the law of large numbers and central 
limit theorem under sublinear expectations, which indicate that G-normal dis- 
tribution plays the same important role in the theory of sublinear expectations 
as normal distribution in the linear expectations (see [12] and [T3]). 

Since the importance of G-normal distribution and G-Brownian motion, 
Peng proposed the problem of how to calculate the G-expectation of <j){B t )l 
For convex or concave </>, Peng gave the formula in [10] and [11] . But for neither 
convex nor concave <j>, how to calculate? In particular, the calculation of B^ n+ , 
for each integer n > 1. 
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In this paper, we give the relation between the solution of G-heat equation 
{T]) with the initial condition <fi(x) — x n for each integer n > 1 and the solution 
of ordinary differential equation (see Section 3), then we can get the solution of 
G-heat equation ([T]) by solving ordinary differential equation. In particular, we 
get the G-expectation of B^ n+1 for each integer n > 1. In fact, we also get the 
solution of the Barenblatt equation with the same initial condition. Finally, we 
point out the application of our result in mathematical finance. 

This paper is organized as follows: in Section 2, we recall briefly the notions 
of G-normal distribution, G-expectation and G-Brownian motion. Also, we 
prove that the G-expectation of B\ n+X is strictly bigger than which is different 
from linear case. The main theorem in which we get the solution of G-heat 
equation (TT]) with the initial condition (j>(x) — x n for each integer n > 1 is 
stated and proved in Section 3. In Section 4, we consider the application of our 
result in mathematical finance. 

2 G-Brownian Motions under G- expect at ions 

Let Q, = Go(M + ) be the space of all K-valued continuous paths (u>t)tes.+ ) with 
u>q = 0, equipped with the distance 

oo 

p(u\u 3 ) := ^2- i [(max |^ - <4\) A 1]. 

i— 1 

The corresponding canonical process B t (oj) — oj t , t 6 [0,oo), for uo € ft. We 
denote G/.Li P (K") the linear space of (local Lipschitz) functions cf> satisfying 

\<Kx)-<Kv)\ <C(l + \xr + \yn\x-yl Vz,yer\ 
for some G > 0, m £ N depending on <f>. 

For each fixed T € [0, oo), we set 

L%(?t) :={<l>{B tl) >>- ,B t J :VneN, t u --- ,t n G [0, T], V</> € C LLip {R n )}. 

It is clear that L° lp (T t ) C L° lp (T T ), for t < T. We also set 

OO 

n=l 

Following [TU] and [TT], we can construct a consistent sublinear expectation 
called G-expectation E[«] : L^T) >— > R satisfying, for each X,Y £ L? p (J r ), 

(a) Monotonicity: if X > Y, then E[X] > E[Y}. 

(b) Constant preserving: E[c] = c, Vc € R. 

(c) Sub-additivity: E[X + Y] <E[JC]+E[F]. 
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(d) Positive homogeneity: E[XX] = XE[X], VA > 0. 

Under G-expectation E[-], the canonical process {B t : t > 0} is called G- 
Brownian motion, and the distribution of B\ is said to be the G-normal distri- 
bution, i.e., for each <ft G Ci.Lipffi), the function 

u(t, x) := E[<j>(x + VtBi)], (t, x) G [0, oo) x R 

is the viscosity solution (see [3]) of G-heat equation (JTJ). Moreover, G-Brownian 
motion has independent and stationary increments, i.e., for each n G N, 4> € 
G; j Li P (R"), s, t > and ii, . . . , t n -\ G [0, i], we have 

E[0(B tl , . . . , Bt„_ 1 , B t+S - B t )] = ±[<p{B tl , . . . , B tn _J] , (2) 

where . . . , x n -i) = E[4>(x\, . . . , x n —i, y/sBi)]. Specially, for each <f> G 

G ; , Llp (R), E[cf)(B t )} = E[^(y/tBi)], yt>0. In particular, for each integer n > 1, 
E[P t 2n+1 ] = t n +*i[fl2n-H]. 

Remark 1 For a G (0,1], G-heat equation (QP is a uniformly parabolic PDE 
and G is a convex function, then it has a unique G 1,2 solution (see pV and fl&j)- 

In the following, P 1 ^ denotes Wiener probability measure on f2, 1?^ always 
denotes the linear expectation with respect to P w . Under P , the canonical 
process {B t : i > 0} is the classical standard Brownian motion. 

Lemma 2 For each fixed a G [0, 1), we have 

(i) For each <\> G C l<Lip (R), E[<f>{B t )} > sup^^ £7 w [0(i/B t )]. 

(ii) For each integer n>l, E[P 2n+1 ] > 0. 

Proof. Noting that u(t,x) := E[(j>(x + B t )} and u"(t,x) := E w [<f>(x+uB t )] are 
respective the viscosity solution of G-heat equation (TT]) and u" — \v 2 u u xx = 
with the same initial condition, by comparison theorem for parabolic partial 
differential equations (see [3]), we get (i). It follows from ([2]) that 

E[B 2n+1 ] = E[(Bi + B x - Bi ) 2n+1 ] = E[0(Pi )], 
where 4>{x) = E[(x + Bi) 2n+1 ]. As a consequence of (i), 
<f>(x) > sup E w [{x + uBi_) 2n+1 ] 

a<v<\ 2 

> E ( 2 " 2 + % 2{n -' l)+1 E W [Bf] + n(2n + 1) (l - a 2 )(x-) 2n ~\ 
t=o V 1 / 2 

Therefore, 

E[B 2n+1 ] = E[0(Bi )] > ^[0(^1)] > n(2n + 1) (l - a 2 )E w [{B- 1 ) 2n - 1 ]. 

2 2 2 2 

As a < 1, we get E[B 2n+1 ] > for each integer n > 1. □ 
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Remark 3 For convex (resp. concave) (f> £ Ci,Lip(JSL) } we have E,[(j)(B t )] — 
E w [cj)(B t )] (resp. t[4>{B t )] = E w [<j){aB t )]); for odd function <j) G C lMp (R), we 
have E[<j)(Bt)] = E[-0(S t )], all see Proposition 7 in JTQ}. For a < 1, t > 0, 
E[_B t 2ll+1 ] = E[— Bf n+1 ] > 0, this is different from linear case, and at the same 
time show that the inequality in (i) can't be changed into equality. 

Remark 4 We define T t := a{B u ,0 < u < t} V N w for each t>0, where N w 
is all P w -null sets, V is the set of all Tt-adapted controls with values in [a, 1]. 
Using a stochastic control method (see ^j), we have for each <f> € Ci t Lip(M), 

= sup£^( / v s dB s )]. 
v.ev Jo 



3 Main Result 

In this Section, we discuss the solution of G-heat equation ([I]) with the initial 
condition cj>(x) — x n for each integer n > 1. Noting that for convex (resp. 
concave) cf>, u(t,x) := E w [<f>(x + B t )] (resp. u(t,x) := E w [<f>(x + aB t )]) is the 
solution of G-heat equation ([T]). Consequently, we only consider the following 
G-heat equation: 

where a e [0, 1], (t, x) <G [0, oo) x R. For a = 1, G-heat equation is classical 
heat equation, then u(t,x) := E w [(x + B t ) 2n+1 ] is the unique solution. In order 
to give the solution of G-heat equation ^ with a S [0, 1), we define for each 
integer n > 0, 



q( x)-Y (2n + 1)! x^ 



(n + i)\(n - i)\ ( (2n + l\ (2n + l 



i=0 



2{n - i))\\{2i)\ \\ / \n-i 



(4) 



„2» 



with the convention that 0! = 1 and 0!! = 1. Now, we give our main Theorem. 

Theorem 5 For each integer n > 1, <?«(•) and /i„(-) defined in The fol- 

lowing statements hold. 
(i) For each fixed a £ (0, 1), we define 



9n(x) + j^[hn(x)exp(-* r )-g n (x)J x exp(-\)dt], x > c n 

o- 2n+1 9n(%) + ^frrMf ) ex P(~^) + 5n(f ) /JL exp(-4) di ], x < c 

(5) 
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where c n , k n and d n are constants such that 



h n -i(—)+g n -x(—)exp(-^) / exp(- — )dt 
a a 2a z J_ 00 2 

c 2 f°° t 2 
= o 2n [h n ^i(c n ) -g n -i(c„)exp(-|^) / exp(-— )dt] 



(2n)!! ffn _ 1 (c„) (6) 



/i n -i(c n )exp(-^-) - g n -i(cn) exp(—^)dt 

rfn = - ^— fc n exp ( — = c n ) 

.9«-i(c„) 2o £ 

k c„ < 

i x 

TTien x) := P°( — ■=) is the unique solution of G -heat equation 

yt 

(ii) for cr = 0, we define 

{k n x 2 f°° t 2 

9n{x) + j^^[K{x)ex.Y>{-^-)~ g n (x) J exp(-—)dt],x>c n 
X ^ X < ^ 

(7) 

where c n , k n are constants such that 

( 2 r°o ,2 

(2n- 1)! = c 2 "[^-i(c„) - .g„_i(c n )exp(^) j exp(-~)di] 

fc« = - ( 2n 2 _ i)n g "- l(c " )c "" ex P(y ) ' S 

C n < 

XTien u n (t, := t n+ ^ P n {—j= ) is the unique viscosity solution of G-heat equation 

yt 

W- 

In order to prove the main Theorem, we need the following Lemmas. The 
first Lemma gives the relation between the solution of G-heat equation ([3|) and 
the solution of the following ordinary differential equation (ODE for short): 

{y")+ - o-VT + xy' - (2n + l)y = 0, lim t n+ ?y{^=) = x 2n+1 . (9) 

t^0+ yt 

Lemma 6 Let a G (0, 1] be a fixed constant. The following statements hold. 

(i) If u is the solution of G-heat equation then P{x) := u(l,x) is the C 2 

solution of ODE 

(ii) If P is the C 2 solution of ODE then u(t,x) := t n+ ^P(^) is the 

yt 

solution of G-heat equation |3J). 
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Proof. For (i), noting that B\ is G-normal distributed, then for t > 0, 

u(t, x) = E[(x + VtB^ 1 } = t n+ ^t[{^p + B^ 1 } = i n+ = u(l, -^=). 

v i v t 

Consequently, u(t,x) = t n+ ^P(—=) for t > 0. Also u(0,x) = £ 2 ™ +1 , we get 

Vi 

hm t _ >0 + i" + 2P(-^) = x 2n+1 . According to Remark 1, we have P e C 2 . It is 
easy to verify that, for t > 0, 

| = („ + i,P(i)-^P'(i)], £ = *»-4P»(i). (10) 

Substituting (I10|) into G-heat equation ©, we conclude the result. The proof 
for (ii) is similar. □ 

Remark 7 Let u a denote the viscosity solution of G-heat equation (0) with 
a £ [0, 1], then for each fixed (t, x) £ (0, oo) x R, </>((t) := u tT (t, x), a 6 [0, 1], zs 
a continuous and decreasing function (see and JlOjj). 

In the following, we first solve ODE © for a £ (0, 1], then by Lemma (6) we 
get the corresponding solution of G-heat equation By the above Remark, 
taking a \ 0, we get the solution of G-heat equation ([3]) with a — 0. For this, 
we give the following Lemmas. 

Lemma 8 Let a ^ 0, a > 0, if ip(-) is a solution of ODE y" + xy' — ay = 0, 
then tp(x) :— <p(~) is a solution of ODE a 2 y" + xy 1 — ay = 0. 

Proof. It is easy to verify the result. □ 

In the following, we use g n (-) and h n (-) defined in 

Lemma 9 The general solution of ODE y" + xy' — (2n + l)y = is 

x 2 f°° t 2 

y{x) = \ig n {x) + \2[h n (x)exp( — — ) - g n (x) J exp(-— )dt], 

where Ai and A2 are arbitrary constants. 

Proof. Applying Lemma 6 with a = 1, we get y n (x) = E w [(x + Bi) 2n+1 ] and 
Vn(x) — V2nE w [((x + i?i)~) 2ll+1 ] are two linear independent solutions of ODE 
y" + xy' — (2n + l)y = 0. It is easy to check that y n (x) = g n (x) and 

x 2 f°° t 2 

y n (x) = K{x) exp( — — ) - g n (x) J exp(-— ) dt, 

which completes the proof. □ 
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Lemma 10 Let x > 0, then for each integer n > 1, we /iai>e 

Proof. We define m n {x) = V2nE w [((x + £i)-) 2 ™ +1 ], then by the proof of the 
above Lemma, we know 

x 2 f°° t 2 

m n {x) = h n (x)exp( — — ) - g n (x) J exp(- — )dt. 

Noting m n (x) > 0, the right-hand side inequality of (JTTJ) holds. We also observe 
that m n (-) is a decreasing function, then m' n (x) < 0, which yields the left-hand 
side inequality of (fTTj) . □ 

Remark 11 It is easy to verify that, for each fixed x > 0, asn-> oo, 
K-l( x ) +9n-i(x) 



„2 . /•oo ,2 



ar I / r 
cx P(-y)| / cxp(-— )<ft. 



Furthermore, we can also prove that 

[°° , t\.. CiM+gn-iW r • r v^(2fc)!(n~fc)! 1 

i. 6XP( ~ 2 +x gB - 1 (x) CXP( "T ) < ^23M } ' 

hJx) . x 2 , f°° . t 2 , , . rV ^7r(2Jfc + l)!(2(n-ife)-l)!n!, 

-^f^exp — -)- / expf- — )dt< min { . W^rrr w r r^l 

5 „(x) vy 2' f x vy 2' i<k<n-i l 2 k x 2 ^ k+1 ){2n)\{n-k-l)\ 1 

Lemma 12 Let a G (0, 1) be a fixed constant, and we define for n > 1, 

X X x 2 f ® t 2 

f n (x) :=h n -i(—) +3„_i(-)exp(— -5-) / exp(- — ) dt - a 2n h n ^ x {x) 
a a 2a l J_ OQ 2 

x 2 f°° t 2 
+ CT 2n g„-i(a;)exp(y) J exp(- — ) dt, 

then there exists a unique Xq < such that f n {%o) = 0- 
Proof. It is easy to verify that 

/ n (0) = (1 - a 2n )(2(n - 1))!! > 0, /„(-oo) = -00, 

2 — 2 

fL(x)=-(9' n -i(-) + -9n-i(-))eM^) \" exp(~)dt 
a a a a iu z J _ oc 2 

+ -(^-l(-) + 9n-l{-)) - <7 2 "(ClW+5n-lW) 

a a a 

, x 2 f°° t 2 

+ o 2n Wn-i{ x ) + x 9n-i{x))e-xv{—) J exp(-—)dt. 
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As x < 0, 

2 — 2 

fn( x ) >-(5n-i(-) + -5n-i(-))exp(^2) / exp(--)df 

+ -(Ki-l(-)+9n-l(-))- 

a a a 

Applying Lemma 10, we get f' n (x) > for x < 0, which completes the proof. □ 

Lemma 13 For each integer n > 1, we have 

x 2 f x t 2 
lim x 2n [h n ^{x) + 5n _i(a:)exp(— ) / exp(--) dt] = (2n - 1)!. 

Proof. It is easy to prove that 

h n -i(x)g n (x) - g n -i(x)h n (x) = (2n - l)lx, 

h n -x(x)g' n _ 1 (x) + xh n ^i(x)g n ~i(x) - h' n _ 1 (x)g n -i(x) - g^-xix) = (2n - 1)!, 
then for x < 0, by Lemma 10, we obtain 



2 



x n [h n -x(x) + g n -i(x) exp(— ) / exp(-— )cft]> 



t 2 N , (2n- l)!a: 2n+1 



2 ,•■< ,2 (2n-l)!ir 2 ™ 



r t 2 ( 

x 2n [h n ^(x) + 3 „_ 1 (a;)exp(— ) / exp(--) dt] < 



(x) + xg n -x(x) ' 



Taking a; — * — oo yields the result. □ 

Proof of Theorem 5. For each fixed a £ (0, 1), we denote by P%(-) the C 2 
solution of ODE ©. We assume that there exists a constant c n < such that 

(tWXz) > for x > c„; (A-P^)(a) < for z < c„. 
do; da; 

Under this assumption, by Lemma 8, Lemma 9 and lim t ^ + t n+ ^ P%(-^) — 
x 2n+1 7 it follows that P^(-) has an expression as in ((5|), where c„, k n and d n are 
undetermined constants. Note that P% € C 2 , we have 

lim P-(x) = lim PI (x), lim ( A_P„ ff )(x) = 0, lim (-^P*)(*) = 0. 

x — >Cn x — 'C n x — >Cn ^X x — >c n ™X 

This yields ([6]) by direct verification. Applying Lemma 12, we know that there 
exists a unique c„ < satisfying the first equation in |6]), then k n and d n are 
also unique. It is easy to verify that P° determined by (0 and ([6]) belongs to 
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C 2 . Now, we show that this P% satisfies the above assumption. For x > c„, 



1 



(2n+l)(2n) y dx 2 ~ 

k n x 2 f°° t 2 

= 9n-i(x) + J2^( h n-i( x ) CX P( — 2") _ 9n-x{x) J exp(-— ) dt) 

_ plfff. , r \2n-li gtCCg-Kgl) ) 2 " ] , n \2n-ll 

~ E [{x + Bl) ] ~ ww^WW^] i^ + B,) } 

— wW\((~ , d \+\2n-li E W [((x + Bi) ) 2n X ] „W\((„ , r \+\2n-ll 

-25 [((z + iJi) ) ] - EwJ^^y)^] E [« C » + Bl ) ) I" 

It is easy to see that, for x > c„, 

B w [((c + Si)-) 2 "- 1 ] > £ w [((z + Si)-) 2 "" 1 ], 
E w [((x + B 1 ) + ) 2n - 1 ] > E w [{{c n + Si)+) 2 "- 1 ], 

therefore, we get (4-zP%)(x) > for x > c n . Similarly, we can prove that 

(■& P n)( x ) ^ for x ^ c «- Tllus > this P « is indeed the C 2 solution of ODE 
((9]). By Lemma 6, u„(t,x) := t n+ ?P%(-^) is the solution of G-heat equation 
ijni, and then the proof of part (i) is complete. We now prove part (ii). We 
define 



l n (x 

It is easy to prove that 



/ exp(- — )dt — exp(- — ), x<0. 

Jx 2 9n-i\x) 2 



(2n- 1)! 



x^O 



lim l n (x) = +00; l' n (x) = ^' ^ exp(-y) > 0, x < 0. ( 12) 



According to Remark 7, u ff (l,0) = k n {a) = L j 2 "j^ j as a T 1, then by fTS) 

we get c„(cr) f as cr j 1. Hence. 
13, we obtain (ii). □ 



{ —00 as cr 4 0. Taking cr J, 0, by Lemma 



Remark 14 For eac/i /lied cr e (0,1), E[ J B 2 ™ +1 ] = <(i,0) = £;„i n+ ^ where 
k n is determined by (iSj). In order to get k n , we solve the first equation in {6]) 
of c n , then k n is uniquely determined by c n . Specially, ~E[Bf] =kt^, where k is 
determined by the following equations: 



1+ r xp( ^ ; 



t 2 

cxp(-— )cft 

1 * 



= cr 2 [l-cexp(y) 
2c 



exp(-Tj-) - c J c °° exp(-^-) 



c < 



cx P(-y)*] 



(13) 



9 



For a — 0, E[B 2n+1 ] = u n (t,0) = k n t n+ ? , where k n is determined by @J. In 
particular, E[I?f] = ktz , where k satisfies the following equations: 



3 r c \ r 

c exp(— ) J exp 



-2c 3 exp( — ) 



c < 



t 2 



(14) 



In /act, for each fixed a G [0, 1], we can get 
corresponding solution. 



IB 2 ] for each I G K by the 



Remark 15 Our method can be extended to the functions 

cj>{x) = h(x+r + i 2 ( x -) a , 

where a > is any given number, l±, l 2 G M. In fact, this is the only type 
of functions satisfying <p(\x) = \ a (j)(x), VA > 0, x G K. It is easy to obtain 
the related Lemma 6, then we get a method to solve G-heat equation (J^l with 
this kind of initial conditions. Specially, for each a G [0, 1], let u° n denote the 
solution of G-heat equation then u°(t,x) := u°(t,—x) is the solution of 
G-heat equation (QP with the initial condition 4>{x) = —x 2n+1 . 



4 Application to Mathematical Finance 

Let {Q,,T,P) be a probability space and (B t )t>o be a 1-dimensional standard 
Brownian motion in this space. We denote by (Ft)t>o the natural filtration 
generated by B, i.e., 

T t := a{B s : s < t} V W, 

where AT denotes the set of all P-null subsets in T . For given a G [0, 1], we 
denote by U a the space of all ^-adapted controls with values in [a, 1]. For fixed 
T > 0, let (St)t<T be the price of a stock and let (log S t )t<T satisfy the following 
stochastic differential equation 

d log S t = \idt + a t dB t , S = 1, 

where /i is a constant, cr. G U a . In mathematical finance, we often need to 
calculate (see [5]) 

sup E P [<p(S T )]. 
<?.eu a 

Specially, for (p(Sr) = (log St)", n > 1, according to Remark 4, 

sup Bp[(logS r )"]=E[(/iT + J B T )"]. 
<y_eu„ 

Therefore, by Theorem 5, we can get sup^ eUa Ep [(log St)" - ] • 
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